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ABSTRACT 


Stress,  particle  velocity,  and  acceleration  data  were 
obtained  from  a  140-ton  nitromethane  row-charge  detonation 
in  a  saturated  clay-shale.  The  charge  row  consisted  of  two 
40-  and  three  20- ton  spheres  (first  and  fourth  charges  were 
40  tons)  extending  S  11  W  from  the  Pre-GONDOLA  I  Charlie 
crater  on  an  80-foot  spacing  (first  charge  105  feet  from 
crater  center) .  The  primary  instrument  line  (stress  and 
particle  velocity  gages)  extended  S  79  E  from  the  center 
charge  (G)  with  gage  stations  100,  150,  200,  300,  and  400  feet 
from  G.  Gage  depths  were  10  and  50  feet  except  the  400-foot 
station  which  was  gaged  at  50  feet  only.  Acceleration  measure¬ 
ments  were  made  on  the  slopes  of  the  Pre-GONDOLA  I  Delta  and 
Charlie  craters  and  at  the  Control  Point  (CP) . 

Peak  stress  measurements  were  higher  than  Pre-GONDOLA  I 
data.  As  the  charge-to-gage  distance  increased,  the  individual 
compressional  waves  from  the  various  charges  began  to  coalesce. 
The  resulting  stresses  (coalesced  wave)  scale  using  combined 
charge  weights  and  cube  root  scaling. 

Initial  peak  horizontal  velocity  compares  well  with 
Pre-GONDOLA  I  data.  Similar  to  stress,  velocities  resulting 
from  the  coalesced  shock  front  scale  using  cube  root. 
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A  low- amplitude,  long-duration  acceleration  was  recorded 
at  the  Delta  and  Charlie  craters.  The  fractured  crater  material 
filtered  out  the  high-frequency,  high- amplitude  motions 
anticipated  at  the  crater  slope  and  lowered  the  average  prop¬ 
agation  velocities.  CP  accelerations  compare  well  with  predic¬ 
tions  based  on  competent  media. 


PREFACE 


This  study  was  performed  by  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  for  the  U.  S.  Army  Engineer 
Nuclear  Cratering  Group  (NCG)  during  the  period  May-June  1967. 
Field  support  was  provided  by  the  Fort  Peck  Area  Office,  CE, 
and  the  Omaha  District,  CE. 

The  work  was  conducted  by  personnel  of  the  WES  Nuclear 
Weapons  Effects  Division  (NWED)  under  the  direction  of 
Mr.  G.  L.  Arbuthnot,  Jr.,  Division  Chief;  Mr.  L„  F,  Ingram, 
Chief,  Physical  Sciences  Branch;  and  Mr.  J.  D„  Day,  Chief, 

Blast  and  Shock  Section.  NWED  Project  personnel  were  Messrs.  C 
Joachim  and  W.  M.  Gay.  Messrs.  L.  T.  Watson  and  G.  H.  Williams 
of  the  Instrumentation  Branch,  Technical  Services  Division,  WES 
were  responsible  for  the  instrumentation.  COL  John  R.  Oswalt, 
Jr.,  CE,  was  Director  of  the  WES  and  Mr.  J.  B.  Tiffany  was 
Technical  Director  during  this  work. 
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CHAPTER  1 


INTRODUCTION 

1.1  DESCRIPTION  OF  PROJECT  PRE-GONDOLA  II 

Project  Pre-GONDOLA  II  was  a  row- charge  cratering  experi¬ 
ment  in  weak,  wet  clay-shale  conducted  by  the  U.  S.  Army 
Engineer  Nuclear  Cratering  Group  (NCG)  as  a  part  of  the  joint 
Atomic  Energy  Commission-Corps  of  Engineers  nuclear  excavation 
research  program.  The  primary  purpose  of  this  nominal  140-ton 
5- charge  row- charge  experiment  was  to  gain  row- charge  cratering 
experience  in  a  weak,  wet  medium.  In  addition,  this  experiment 
tested  techniques  for  connecting  a  row- charge  crater  to 
an  existing  crater  and  for  over-excavating  to  accept  throwout 
from  a  follow-on  connecting  row- charge  crater. 

Project  Pre-GONDOLA  II  was  detonated  in  Valley  County, 
near  the  edge  of  the  Fort  Peck  Reservoir  approximately  18 
miles  south  of  Glasgow,  Montana,  at  exactly  0800  hours  (MDT) 

28  June  1967.  Coordinates  of  the  center  charge  were  W106  38' 3l" 
N47  55'5l''.  The  orientation  of  the  row  was  along  a  10  east 
of  north  to  10  west  of  south  alignment  extending  through 
the  center  of  the  Charlie  crater. 

The  average  lip- crest  to  lip-crest  dimensions  were  640 
feet  by  280  feet.  Individual  charge  yields,  depths,  spacing, 


and  resulting  apparent  crater  dimensions  were  as  follows: 

Apparent 


Emplacement  Configuration 

Crater  Dimensions 

(feet) 

Depth  of 

(feet) 

Charge 

tons  NM** 

Burst 

Spacing 

Width 

W  Depth  D 
—  a  - a 

Charlie  (existing 
crater) 

105.5 

E  (connecting 

38.61 

59.7 

206.5 

55.5 

crater) 

79.8 

F 

19.70 

49.4 

79.9 

152.5 

37.5 

G  (center  crater) 

19.55 

48.8 

79.9 

164.0 

36.9 

H  (over-exca¬ 

39.56 

59.9 

214.5 

57.0 

vating  charge) 

79.9 

I 

20.00 

48.8 

173.0 

33.5 

1.2  OBJECTIVES 

The  primary  objective  of  the  study  described  herein  was 
to  obtain  close-in  (<400  feet)  ground-motion  and  earth- stress 
measurements  from  a  row- charge  detonation  in  weak,  wet  clay- 
shale.  Motion  measurements  were  also  performed  for  the  existing 
Charlie  and  Delta  crater  slopes  and  at  the  Control  Point  (CP) „ 


1 . 3  BACKGROUND 

Efficient  use  of  nuclear  devices  for  large-scale  exca¬ 
vations  requires  reasonably  precise  predictions  of  initial 

**  Spherical  charges  of  liquid  explosive  nitromethane  (CHgNO^ . 
Actual  total  charge  weight  was  137.42  tons. 
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crater  dimensions  and  subsequent  long-term  stability.  Such 
predictions  can  only  be  made  if  cratering  experience  is  gained 
in  a  wide  variety  of  geologic  media,  particularly  those  which 
might  be  encountered  in  excavation  projects.  The  physical 
measurement  of  crater  parameters  together  with  measurements 
of  transient,  directly  coupled  ground  motions  and  stresses 
is  used  to  evaluate  the  engineering  properties  of  the  excavated 
media . 

Early  cratering  experiments  were  conducted  by  NCG  in  dry 
materials,  i.e.,  basalt,  rhyolite,  and  desert  alluvium.  Since 
it  is  highly  probable  that  some  portions  of  any  major  excavation 
will  be  in  material  at  least  partially  saturated,  cratering 
data  are  needed  for  wet  materials.  Also,  design  for  long¬ 
term  use  of  a  canal  requires  knowledge  of  slope  stability 
when  explosively  excavated  slopes  are  subjected  to  wetting. 

This  includes  pore  pressure  effects  of  water  behind  the 
slopes  and  erosive  effects  of  water  flowing  in  the  channel. 

The  Pre-GONDOLA  cratering  experiments  conducted  by  NCG  at 
Fort  Peck,  Montana,  provide  data  for  one  such  medium,  a  weak, 
wet  clay- shale  (known  as  "Bearpaw"  shale) . 

Due  to  a  lack  of  meaningful  cratering  data  for  these  media, 
it  was  necessary  to  conduct  a  site  calibration  test  series. 
Pre-GONDOLA  I,  a  series  of  four,  single-charge  20-ton  experiments, 
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was  conducted  for  this  purpose.  Close-in  stress  and  ground 

motion  measurements  were  made  in  conjunction  with  the  Bravo 

1* 

Event  of  this  series  .  The  data  obtained  served  as  a  basis  for 
design  of  this  project  on  Pre-GONDOLA.  II,  a  row-charge  event. 
This  report  describes  that  portion  of  the  Pre-GONDOLA.  II 
experiment  concerned  with  measurement  of  close-in  stress  and 
ground  motion. 


*  Superscript  numerals  refer  to  similarly  numbered  items 
in  References  at  the  end  of  the  text. 
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CHAPTER  2 


PROCEDURE 

2.1  EXPERIMENTAL  PLAN 

The  main  gage  layout  for  Pre-GONDOLA  II  is  shown  in 
plan  and  profile  in  Figure  2.1.  A  total  of  nine  earth 
stress  gages,  and  nine  horizontal  and  seven  vertical  particle 
velocity  gages  were  installed  at  nominal  depths  of  10  and  50 
feet  and  at  horizontal  distances  of  100  to  400  feet.  The  deep 
gages  (50-foot  depth)  were  located  along  a  line  perpendicular 
to  the  charge  row  at  charge  G.  The  shallow  gages  (10-foot  depth) 
were  installed  in  separate  holes  located  10  feet  north  of 
and  perpendicular  to  the  main  gage  line.  In  addition  to  the 
main  gage  line,  surface  accelerometers  were  placed  on  the 
existing  Charlie  and  Delta  crater  slopes  (Figure  2.2)  to 
measure  motion  of  the  slopes,  and  surface  accelerometers  were 
placed  at  the  CP.  Table  2.1  lists  the  main  line  instrumentation 
by  gage  type  and  location  and  also  gives  predicted  peak  values. 
Table  2.2  lists  all  accelerometers,  and  gives  peak  predicted 
values . 

2.2  PREDICTIONS 

Because  of  the  uniqueness  of  the  experiment  there  was  no 
empirical  basis  for  predicting  stresses  and  motions.  Although 
close-in  earth  motions  from  multiple  contained  bursts  in  a 
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linear  array  have  been  measured  ,  these  measurements  were  made 
in  an  entirely  different  material  (dry  basalt)  and  in  line 
with  the  charges  rather  than  normal  to  the  array  as  planned 
for  Pre-GONDOLA  II.  Therefore,  predictions  were  based  on 
stress  and  velocity  measurements  obtained  from  the  20-ton 
single-charge  calibration  series  (Event  Bravo).  Based 
on  these  measurements,  stress  and  horizontal  particle  velocity 
inputs  from  the  individual  charges  of  the  row  were  computed 
for  each  of  the  nine  gage  locations.  Cube- root  scaling  from 
Bravo  was  used  to  obtain  the  40-ton  charge  inputs.  Applying 
the  principle  of  superposition,  peak  values  which  occurred 
within  any  10-msec  time  interval  were  added.  The  maximum 
values  thus  obtained  were  used  for  instrument  set-range 
predictions  (Table  2.1).  These  predictions  were  considered 
conservative  because,  in  this  procedure,  no  allowance  or  con¬ 
sideration  was  made  for  the  signal  decay. 

Peak  particle  velocity  predictions  were  based  on  radial 
measurements  from  Pre-GONDOLA  I1.  Using  these  data,  vertical 
and  horizontal  particle  velocity  predictions  were  computed 
from  the  following  relations: 

u  =  u  cos  0 
H  R 

■  UV  =  UR  Sin  6 
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where 


u 

u 

u 


R 

H 

V 


e 


=  radial  particle  velocity 
=  horizontal  particle  velocity 
=  vertical  particle  velocity 

=  angle  between  horizontal  and  charge-gage  line 


Particle  velocities  up  to  1/3  u  have  been  measured 

K 

perpendicular  to  the  compressional  wave  propagation  although 
wave  propagation  theory  does  not  predict  these  motions.  Therefore, 
the  minimum  vertical  particle  velocity  prediction  used  was 


1/3  V 

No  acceleration  data  were  available  on  crater  slopes; 

therefore,  predictions  (Table  2.2)  were  made  using  the 

3 

following  equation  taken  from  Newmark  ,  based  on  continuous 
media . 


where 


A 

u 

t 


R 

R 

R 


=  radial  acceleration 
=  radial  particle  velocity 

=  rise  time  to  peak  radial  particle  velocity 


This  equation  is  derived  for  rock,  primarily  from  tests  conducted 
at  the  Nevada  Test  Site. 

The  peak  radial  particle  velocity  obtained  from  the  20-ton 
Bravo  Event  was  used  in  this  equation.  It  was  assumed  that  the 
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single  charge  nearest  the  gage  would  contribute  the  maximum 
ground  motion  input  at  that  gage  location  due  to  the 


unsymmetrical  gage-charge  geometry  (Figure  2.2). 

The  range  of  rise  times  given  by  Newmark  is  1/6  to  1/12 
of  the  arrival  time.  A  rise  time  of  one- twelfth  the  arrival 
time  was  used  in  making  these  predictions.  This  method  gave 
the  highest  peak  acceleration  prediction  and  was  the  more 
conservative  approach.  Peak  vertical  accelerations  were  predicted 
using  the  same  method  used  for  peak  vertical  particle  velocity. 

2.3  INSTRUMENTATION 

2.3.1  Gages.  Earth  stress  measurements  were  made  with 

1/4-  and  1/2- inch-diameter  tourmaline  gages  purchased  from 

Crystal  Research,  Inc.,  Cambridge,  Massachusetts.  The  1/2-inch 

gages  were  employed  in  the  predicted  stress  region  from 

2,500  psi  to  3,500  psi,  and  the  1/4-inch  gages  for  pressures 

higher  than  3,500  psi.  These  gages  had  been  used  successfully 

.  1 

in  previous  tests  in  the  Pre-GONDOLA  series  . 

Velocity  gages  used  in  this  study  were  Sandia  Model  DX  gages 
purchased  from  Spartan  Industries,  Albuquerque,  New  Mexico. 

These  gages  operate  on  the  principle  of  the  overdamped  acceler¬ 
ometer.  A  variable  reluctance  sensor  is  used.  A  more  detailed 
gage  description  can  be  found  in  Reference  1.  Vertical  particle 
velocity  gages  used  were  horizontal  gages  modified  by  the 
addition  of  a  spring  to  counteract  gravitational  forces  on  the 


pendulum . 
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Accelerometers  used  were  of  two  types:  strain  gage 
and  variable  reluctance.  The  accelerometer  designated  "200  AR  CL," 
was  a  1 , 000-G- capacity  instrument  made  by  Consolidated 
Electrodynamics  Corporation,  Pasadena,  California,  and 
that  designated  "150  AR  CS”  was  made  by  Endevco  Corporation, 
Pasadena,  California.  The  remaining  gages  were  variable 
reluctance  type  made  by  Pace  Engineering  Company,  North 
Hollywood,  California. 

2.3.2  Gage  Installation.  A  horizontal  and  a  vertical 
velocity  gage  were  mounted  in  a  5-1/2- inch- diameter  canister 
with  associated  calibration  electronics.  The  canister 
assembly  was  filled  with  paraffin  for  waterproofing  and 
shockproofing  of  the  contents.  The  electrical  cable  was 
fed  to  the  ground  surface  through  a  3/4- inch-diameter  aluminum 
pipe  threaded  into  the  canister  lid.  The  pipe  also  served 
as  a  placement  and  orientation  device. 

Signal  cable  from  the  tourmaline  stress  gages  was  pro¬ 
tected  with  1-inch  rigid  plastic  pipe  to  the  ground  surface. 

The  gage  protruded  1  inch  below  the  end  of  the  plastic  pipe. 

For  installation  the  stress  gage  assembly  was  taped  to  the 
velocity  gage  canister  so  that  it  extended  6  inches  below 
the  bottom  of  the  canister. 

All  main  gage  line  instrument  holes  were  deeper  than 
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required  gage  depths  and  contained  water  to  within  a  few 
feet  of  the  surface.  These  were  filled  with  sand  to  approxi¬ 
mately  1  foot  below  the  desired  gage  depth.  This  was 
accomplished  by  shoveling  sand  into  the  hole  and  allowing 
it  to  settle  through  the  water.  The  entire  gage  assembly 
was  then  lowered  to  the  desired  depth  and  oriented.  Sufficient 
sand  was  shoveled  into  the  hole  to  ensure  that  the  tourma¬ 
line  pressure  gage  and  part  of  the  canister  were  surrounded 
by  sand.  A  grout  mix  was  then  placed  by  1- inch-diameter 
tremie  pipe  from  the  top  of  the  sand  forming  a  3- foot- long 
plug.  When  the  grout  had  set,  the  remainder  of  the  hole 
was  backfilled  with  sand. 

Horizontal  and  vertical  accelerometers  were  mounted  in 
a  3-1/2-  by  3-1/2-  by  6- inch  aluminum  box  which  was  filled 
with  paraffin.  The  box  was  fastened  to  a  4-  by  4- inch  wood 
post  that  had  been  buried  in  the  crater  wall. 

Electrical  signals  from  the  velocity  gages  and  the 
two  variable  reluctance  accelerometers  located  in  the  Charlie 
crater  were  carried  on  four- conductor  shielded  cable  to  a 
junction  box  (Figure  2.1)  located  approximately  500  feet  from 
the  charge  row.  Connections  between  the  junction  box  and 
the  WES  instrumentation  van  were  made  with  a  50-pair  telephone 
cable.  The  two  strain  gage  accelerometers  required  individual 
shielded  cables  and  were  not  connected  to  the  junction  box. 


14 


Four- conductor  shielded  cable  was  used  for  all  gages  in  the 
Delta  crater. 

Piezoelectric  gage  signals  were  carried  on  coaxial  cable. 

A  piece  of  "low  noise”  cable,  approximately  100  feet  long, 
was  used  at  the  gage.  The  remaining  cable  was  standard  coaxial. 
The  low  noise  cable  was  utilized  to  prevent  cable  noise  from 
obscuring  the  portion  of  the  pressure-time  history  of 
interest.  It  was  not  necessary  to  bury  any  surface  cable  runs 
since  the  data  of  interest  would  have  been  retrieved  before 
damage  from  crater  ejecta. 

2.3.3  Recording  Systems.  All  instrumentation  outputs, 
except  from  the  tourmaline  stress  gages  and  strain  gage  acceler¬ 
ometers,  were  recorded  on  light-beam  galvanometer  oscillographs. 
The  high-frequency  outputs  of  the  tourmaline  stress  gages 
and  strain  gage  type  accelerometers  were  recorded  on  an  FM 
magnetic  tape  recorder.  Signal  conditioning  equipment  varied 
with  the  type  of  instrument;  velocity  gages  and  variable 
reluctance  type  accelerometers  used  a  3  kcps  carrier-demodulator 
system;  strain  gage  type  accelerometers  used  a  solid-state 
DC  operational  amplifier  system;  and  the  tourmaline  stress 
gage  used  a  cathode  follower  system. 

All  electronic  recording  equipment  was  housed  in  the 
WES  recording  van  located  approximately  2,100  feet  south  of 
the  Pre-GONDOLA  II  row  charge. 
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TABLE  2.1  MAIN  GAGE  LINE:  GAGE  TYPE  AND  SET  RANGES 
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Note:  Horizontal  distances  shown  are  to  closest  charge. 


TABLE  2.2  ACCELEROMETER  LOCATION  AND  SET  RANGES 
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Horizontal  distances  shown  are  to  closest  charge. 
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Layout  of  main  gage  line. 
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CHAPTER  3 


RESULTS  AND  DISCUSSION 

3.1  INSTRUMENT  PERFORMANCE 

All  timing  and  sequence- initiating  signals  were  received 
and  properly  translated,  resulting  in  all  recording  equipment 
operating  as  planned. 

Of  the  35  gages  installed  for  Pre-GONDOLA  II,  32  were 
operable  at  shot  time.  Gage  100  S  50  apparently  failed  prior 
to  detonation  although  resistance  readings  had  given  no 
indication  of  this  condition.  Equipment  in  the  WES  van 
sustained  damage  when  the  site  transformer  bank  was  hit  by 
lightning,  and  lightning  damage  is  suspected  in  the  preshot 
failure  of  gage  200  AR  CL,  a  strain  gage  accelerometer.  After 
equipment  repairs,  a  gage  check  revealed  this  gage  to  be 
inoperable.  A  resistance  check,  upon  removal  from  the 
canister,  showed  all  leads  open.  The  cause  of  failure  in 
200  AV  CL  is  not  known.  Three  of  the  remaining  gages 
(100  UV  10,  100  UV  50,  and  200  UV  50)  were  overdriven  and 
rendered  only  qualitative  data. 

3.2  DATA  REDUCTION 

The  oscillograms  and  tape  playbacks  were  read  on  an  electro¬ 
mechanical  digitizer,  and  the  data  converted  to  card  form.  These 
data  cards  vrere  then  processed  through  a  digital  computer 
where  a  temperature  correction  and  an  integration  were  performed 
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on  the  particle  velocity  records  to  yield  displacements.  All 
data  were  then  automatically  plotted. 

Peak  measurements  of  particle  velocity  and  stress  are 
summarized  in  Table  3.1.  Acceleration  peak  values  are  presented 
in  Table  3.2.  Time  histories  of  particle  velocities  and  their 
associated  displacements,  soil  stresses,  and  accelerations 
are  presented  in  Appendix  A. 

3.3  ARRIVAL  TIMES 

Compressional  wave  front  arrival  times  at  all  main  line 
gages  are  plotted  versus  slant  distance  in  Figure  3.1.  The 
average  propagation  velocity  based  on  gages  at  charge  depth 
was  6,600  ft/sec.  Gages  at  a  nominal  10- foot  depth  yielded 
slightly  longer  arrival  times  and,  consequently,  a  slower 
propagation  velocity.  The  higher  velocity  at  depth  is  probably 
due  to  more  competent  material  at  charge  depth  compared  to 
the  weathered  zone  near  the  surface. 

The  average  compressional  wave  velocity  (6,600  ft/sec)  for 
the  main  gage  line  was  lower  than  the  value  reported  for  the 
Pre-GONDOLA  I1  Bravo  Event.  The  Pre-GONDOLA  I  value  was  based 
on  gages  spread  over  a  greater  distance.  The  shock  wave 
traveled  through  deeper,  more  competent  material  to  reach 
the  farther  gage  stations.  The  result  was  a  slightly  higher 
velocity.  The  accelerometers  located  at  the  CP  yielded  a 
velocity  comparable  to  the  Pre-GONDOLA  I  value  of  6,800  ft/sec. 
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More  precise  arrival  times  were  obtained  from  the  piezo¬ 
electric  gage  traces  than  from  the  particle  velocity  records 
because  of  the  ability  of  the  stress  gages  to  respond  faster. 

3.4  SOIL  STRESS 

Initial  peak  stress  versus  slant  distance  is  plotted  in 

1  4,  5 

Figure  3.2.  The  20-ton  Pre-GONDOLA  I  and  scaled  granite 
data  are  shown  in  this  figure  for  comparison.  The  first 

disturbance  sensed  by  all  pressure  gages  in  the  main  gage 

array  was  caused  by  charge  G  (with  the  exception  of  stations 

300  and  400  feet  from  G,  both  at  50-foot  depths).  The 

compressional  waves  from  charges  F,  G,  and  H  coalesced  after 

traveling  approximately  290  feet  and  are  not  separately 

identifiable  300  feet  from  charge  G.  The  compressional  waves 

from  charges  E  and  I  still  lagged  the  coalesced  wave  slightly 

at  the  400-foot  station  (see  Figure  A. 12).  Data  from  the 

deep  gage  stations  300  and  400  feet  from  G  are  presented  at 

both  actual  (80  tons)  and  scaled  distances  (20  tons)  in 

1/3 

Figure  3.2.  The  scaling  factor  is  S  =  (W  /W^)  where: 

=  combined  charge  weight  (F,  G,  and  H) 

W  =  charge  weight  being  scaled  (20  tons) 

2 

Initial  peak  stress  data  from  the  deep  gages  are  consistently 
higher  than  those  from  shallow  gages.  The  deep  gages  were 
located  in  stronger,  saturated,  more  competent  material  which 
transmitted  stronger  shock.  The  shallow  gages  were  apparently 
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located  above  the  groundwater  level  which  is  estimated  to 

be  between  10  and  20  feet  below  ground  surface  in  the  Pre-GONDOLA  II 
6 

area  .  In  addition,  wave  forms  recorded  by  shallow  gages 
at  stations  200  and  300  feet  from  charge  G  (gages  200  S  10 
and  300  S  10,  Figures  A. 6  and  A. 9)  indicate  the  presence 
of  an  unloading  or  tensile  (rarefaction)  wave  which  further 
reduced  peak  values  at  these  stations.  The  gage  station 
150  feet  from  charge  G  yielded  initial  peak  stress  data 
significantly  higher  than  the  general  trend  (gages  150  S  10 
and  150  S  50,  Figure  A. 3).  That  both  deep  and  shallow 
gages  responded  in  this  manner  indicates  some  cause  other 
than  gage  error. 

When  a  compressional  wave  impinges  on  a  free  surface  or 

an  interface  of  lower  acoustical  impedance,  it  is  reflected 

7 

as  a  tensile  or  unloading  wave  .  Tension  is  not  recorded  as 
such  but  is  manifested  as  a  decrease  in  compressive  stress. 

When  a  gage  is  located  in  a  medium  with  a  sufficiently  close 
boundary  such  that  the  difference  between  direct  and  reflected 
ray  paths  is  small,  tensile  relief  resulting  in  lower  peak 
values  can  occur.  The  shallow  gages  were  located  close 
enough  to  the  surface  that  tensile  relief  must  be  investigated. 
Consider  the  following: 
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Gage 


Computed  A  T 
(reflected- direct) 


100  S  10  150  S  10  200  S  10 

(Figure  A .  1)  (Figure  A. 3)  (Figure  A. 6) 

1.5  msec  1.4  msec  0.7  msec 


300  S  10 
(Figure  A  .9) 

0.3  msec 


Apparent  rise  time  0.3  msec  0.9  msec 
(first  peak) 


0.5  msec  0.9  msec 


The  shock  wave  from  a  spherical  source  emanating  into  a 
homogeneous,  continuous  media,  such  as  Bearpaw  clay-shale, 
disperses  and  attenuates7.  This  causes  a  rounding  of  the 
shock  front,  resulting  in  longer  rise  times  with  distance.  The 
comparison  presented  earlier  shows  a  rise  time  increase  with 
distance  out  to  the  150-foot  station  (gage  150  S  10)  while  a 
shorter  apparent  rise  time  was  measured  at  the  200- foot 
station  (gage  200  S  10) .  The  computed  difference  between 
reflected  and  direct  arrival  times  at  the  200-foot  station 
is  very  close  to  the  apparent  rise  time  value,  making  tensile 
relief  at  this  station  a  distinct  possibility.  A  similar 
comparison  at  the  300-foot  station  (gage  300  S  10)  indicates 
that  tensile  relief  had  taken  place,  resulting  in  a  lower 


initial  peak  stress  value. 

Peak  soil  stress  data  are  plotted  versus  slant  distance 
in  Figure  3.3.  One  gage  (100  S  10)  was  damaged  after  recording 
the  initial  peak  stress,  resulting  in  loss  of  additional 
stress  data  from  this  gage.  The  deep  gages  recorded  higher 
peak  values  than  the  shallow  gages.  This  was  partially  due 
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to  the  stronger,  saturated,  more  competent  material  found  at 
depth  but  primarily  was  the  result  of  the  gage- charge  geometry 
used  for  the  shallow  gages.  Placement  of  these  gages  10  feet 
north  of  the  main  gage  line  resulted  in  longer  travel 
times  for  the  stronger  compressional  wave  from  charge  H  (40  tons) 
than  for  charge  F  (20  tons) .  This  difference  decreases  at 
greater  distances.  It  is  particularly  evident  when  comparing 
deep  and  shallow  gage  data  at  the  150-foot  station  (gages 
150  S  10  and  150  S  50,  Figure  A. 3)  as  follows: 

Charge  Gage  150  S  10 _  Gage  150  S  50 _ 

Weight,  Slant  distance,  Travel  time,  Slant  distance,  Travel  time, 
No.  tons  feet _  msec  feet _  msec 


G 

20 

155 

24.1 

150 

21.9 

F 

20 

170 

25.6 

170 

24.4 

H 

40 

182 

26.4 

170 

24.4 

The  wave 

form  at  the 

deep  150-foot  station 

(gage  150  S  50) 

shows  that,  as  expected,  the  charge  F  and  H  compressional  waves 
arrived  almost  simultaneously  but  the  charge  H  compressional 
wave  lagged  the  F  wave  by  almost  1  msec  at  the  shallow  gage 
(150  S  10) .  This  results  in  three  distinct  peaks  (compressional 
waves  from  F,  G,  and  H)  with  only  a  small  increase  in  total 
peak  stress  at  the  shallow  gage. 
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Examination  of  time-histories  from  the  deep  gages  (Figures  A. 3, 

A. 6,  A. 9,  and  A. 12)  indicates  that  the  gaged  area  is  in  a 
transition  zone.  Peak  stress  at  the  150-foot  station  (Figure  A. 3) 
is  the  sum  of  a  partially  decayed  pulse  from  a  20-ton  charge 
(G)  and  the  contribution  of  a  40-ton  (H)  and  a  20-ton  (F) 
charge.  As  the  distance  increased,  these  compressional 
waves  tended  to  coalesce.  The  wave  form  recorded  at  the 
400- foot  station  (gage  400  S  50)  indicates  that  charges 
F,  G,  and  H  have  coalesced  and  input  from  charges  E  and  I 
still  lagged  slightly. 

The  initial  peak  pressure  attenuated  as  the  -2  power  of 
distance  as  was  found  for  the  Pre-GONDOLA  I  series.  Total  peak 
stress  also  attenuated  as  the  -2  power  of  distance  at  the 
300-  and  400-foot  stations.  It  must  be  noted  that  the  effective 
input  continually  changed  over  the  gaged  region  due  to 
coalescence  of  shock  waves  from  the  separate  charges.  Although 
the  compressional  waves  from  all  five  charges  had  not  coalesced 
at  the  400-foot  station  it  is  probable  that  after  coalescence 
attenuation  would  continue  as  the  -2  power  of  distance. 

3.5  PARTICLE  VELOCITY 

Initial  horizontal  particle  velocity  data  are  presented 
in  Figure  3.4.  Here  (as  with  the  stress  measurements)  the 
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initial  peak  values  are  inputs  from  charge  G  (20  tons)  except 
the  deep  gages  (300  UR  50  and  400  UR  50)  at  stations  300  and 
400  feet  from  G.  Data  from  these  gages  are  presented  in  Figure 
3.4  at  both  actual  and  scaled  distances.  An  effective  charge 
weight  of  80  tons  was  assumed  for  scaling.  Wave  forms  recorded 
by  the  shallow  gages  (200  UR  10  and  300  UR  10)  at  stations 
200  and  300  feet  from  G  indicate  that  tensile  relief  has 
taken  place.  In  addition,  the  shallow  gage  at  200  feet  from  G 
was  tilted  with  the  arrival  of  the  initial  compressional  wave 
and  registered  a  higher  than  normal  peak  value.  This  data 
point  is  questionable. 

Pre-GONDOLA  I  data  (horizontal  particle  velocity)  are 
included  in  Figure  3.4  for  comparison.  These  data  were  not 
corrected  for  temperature  in  the  original  report1.  The  gage 
calibration  procedure  was  performed  at  65  F.  Assuming  the 
same  ground  temperature  as  at  Pre-GONDOLA  II  (50  F) ,  a  temperature 
correction  factor  of  1.18  is  applied.  The  three  closest 
Pre-GONDOLA  I  measurements  are  indicated  as  low  in  Reference  1. 

The  gage  located  at  a  slant  distance  of  97  feet  was  overdriven, 
while  those  at  slant  distances  of  100  and  133  feet  experienced 
tensile  relief.  The  remaining  Pre-GONDOLA  I  data  points  compare 
well  within  the  usual  data  scatter  (i.e.  20  percent)  with  the 
Pre-GONDOLA  II  data.  . 

Vertical  particle  velocity  data  are  presented  in  Figure  3.5, 
and  exhibit  considerable  scatter.  Wave  forms  registered  at  the 
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shallow  stations  200  and  300  feet  from  G  (gages  200  UV  10  and 
300  UV  10)  indicate  that  there  has  been  some  tensile  relief  of 
peak  values.  The  gage  at  the  station  200  feet  from  G  (gage 
200  UV  10)  also  indicates  early  tilting,  making  the  initial 
peak  value  questionable. 

Particle  velocity-time  histories  and  first  integrals 
(displacement)  are  presented  in  Appendix  A.  No  base  line 
corrections  were  applied  to  the  velocity  data.  Therefore, 
the  integral-time  history  (displacement),  relative  data  at 
best,  must  be  regarded  as  only  qualitative. 

3.6  ACCELERATION  MEASUREMENTS 

Peak  acceleration  data  are  plotted  versus  horizontal 
distance  in  Figure  3.6.  In  retrospect,  the  peak  acceleration 
predictions  were  overly  conservative.  They  were  based  on 
the  assumption  of  competent  material  and  did  not  take  into 
account  the  presence  of  the  intervening  crater  (in  the  case  of 
measurements  on  the  Charlie  crater  lip) .  Moreover,  only 
empirical  information  was  available  preshot.  The  combination 
of  these  factors  (effect  of  crater  boundaries  and  rubble) 
strongly  influenced  the  maximum  surface  accelerations  at 
the  close-in  stations.  However,  estimates  for  the  CP  motions 
were  much  better,  and  excellent  measurements  were  obtained  at 
this  location. 
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Based  on  arrival  time  data  at  Delta  crater  (gages 
380  AR  DL  and  425  AV  DS) ,  a  propagation  velocity  of  4,000  ft/sec 
was  computed  using  the  distance  to  the  closest  charge  (I) . 

This  value  is  much  slower  than  the  6,600  ft/sec  value  computed 
for  the  main  gage  line  instrumentation.  The  discrepancy  is 
explained  as  follows.  Propagation  velocity  computations  are 
an  average  value  for  all  the  material  between  two  points,  in 
this  case  between  the  closest  charge  (I)  and  the  Delta  crater 
gages.  Two  distinctly  different  materials  are  encountered  by 
a  shock  wave  traversing  this  line;  crater  rubble  and  a  more 
competent  clay- shale.  Consider  the  gage  located  on  the  Delta 
crater  slope  (425  AV  DS)  with  an  arrival  time  of  106.6  msec. 
Assuming  a  propagation  velocity  of  6.6  ft/msec  in  clay-shale 
and  1.5  ft/msec  in  the  rubble,  a  path  approximately  80  feet 
long  through  the  rubble  would  result  in  an  average  velocity 
of  4,000  ft/sec  at  this  gage.  Delta  crater  rupture  zone 
measurements  were  not  made.  Excavation  of  the  Bravo  crater 
(also  a  single-charge  20-ton  experiment  in  similar  material) 
gave  the  following  information*: 


*  PNE-1103,  in  preparation. 
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Radius  from  Charge, 
feet 


Depth  to  Base  of  Rupture 
_ Zone,  feet _ 


125 

60 

175 

43 

225 

21 

375 

0 

Gage  425  AV  DS  was  located  approximately  10  feet  below  the 
original  ground  surface  at  a  50- foot  radius  from  the  center 
of  the  Delta  crater.  Based  on  the  above  information,  an 
80-foot-long  path  in  rubble  appears  realistic. 

The  accelerometer  (150  AV  CS)  located  on  the  Charlie 

crater  slope  was  situated  on  the  far  side  of  the  crater  away 

from  the  closest  charge  (E) .  The  shock  wave  reaching  the 

gage  had  to  travel  under  or  around  the  intervening  crater 

through  the  rubble.  This  resulted  in  a  longer  compressive 

wave  ray  path  to  the  gage  and  a  higher  average  velocity  than 

the  2,000  ft/sec  computed  using  the  horizontal  gage-charge 

distance.  Charlie  crater  is  located  at  the  edge  of  the  Fort 

Peck  Reservoir  which  had  a  pool  elevation  at  Pre-GONDOLA  II 

shot  time  of  2,243*6 .  This  was  5  feet  above  the  pool  elevation 

at  the  end  of  the  Pre-GONDOLA  I  series.  The  bottom  of  the  Charlie 

crater  is  at  elevation  2,220  which  is  18  feet  below  the  pool 

elevation  at  the  end  of  Pre-GONDOLA  I.  Therefore,  the  rubble 
"*  Elevations  are  in  feet  referred  to  mean  sea  level. 
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under  the  Charlie  crater  was  probably  saturated,  which  would 
result  in  higher  average  velocities  at  this  location  than  at 
the  Delta  crater. 
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TABLE  3.1  PEAK  MEASUREMENTS  FROM  VELOCITY  AND  STRESS  GAGES  (MAIN  GAGE  LINE) 
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Note:  Horizontal  distances  shown  are  to  closest  charge. 


TABLE  3.1  PEAK  MEASUREMENTS  FROM  VELOCITY  AND  STRESS  GAGES  (MAIN  GAGE  LINE)  (continued) 
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TABLE  3.2  PEAK  MEASUREMENTS  FROM  ACCELEROMETERS 
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ARRIVAL  TIME,  MSEC 
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Figure  3-1  Arrival  time  versus  distance. 
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PEAK  SOIL  STRESS,  PSI 
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Figure  3.3  Peak  soil  stress  versus  distance. 
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INITIAL  PEAK  HORIZONTAL  PARTICLE  VELOCITY,  FT/SEC 
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INITIAL  PEAK  VERTICAL  PARTICLE  VELOCITY,  FT/SEC 
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Figure  3*5  Initial  peak  vertical  particle  velocity  versus  distance. 
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PEAK  ACCELERATION, 
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Figure  3.6  Peak  acceleration  versus  distance. 
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CHAPTER  4 


CONCLUSIONS 

The  computed  Pre-GONDOLA  II  propagation  velocity  was 
6,600  ft/sec.  This  is  lower  than  the  value  reported  for 
Pre-GONDOLA  I.  The  Pre-GONDOLA  I  value  is  an  average  value 
over  a  longer  distance.  The  shock  wave  reaching  the  farther 
gages  travels  through  a  deeper,  more  saturated,  more  competent, 
denser  material. 

Initial  peak  stress  data  are  the  result  of  input  from 
charge  G  (20  tons)  at  all  stations  except  the  deep  gages  300 
and  400  feet  from  G.  These  exceptions  scale  using  the  combined 
charge  weights  of  F,  G,  and  H  (80  tons).  Wave  forms  indicate 
tensile  cutoff  of  peak  values  at  shallow  stations  200  and 
300  feet  from  G.  Although  higher  initial  stresses  were  measured 
from  Pre-GONDOLA  II  than  I,  both  attenuated  as  the  -2  power  of 
distance . 

Peak  stress  measurements  were  obtained  in  a  transition 
zone  where  the  effective  charge  input  was  changing  with 
distance.  As  the  charge- to- gage  distance  increased,  individual 
compressional  waves  began  to  coalesce,  resulting  in  higher 
stresses  with  distance  than  found  for  initial  peak  data. 

Pulses  from  charges  F,  G,  and  H  coalesced  approximately 
290  feet  from  charge  G,  but  all  five  charge  inputs  had  not 
coalesced  at  the  400- foot  station.  Placement  of  the  shallow 
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gages  10  feet  north  of  the  main  gage  line  resulted  in  a  loss 
of  symmetry.  Thus,  greater  distances  were  required  for  the 
individual  compressional  waves  to  coalesce.  Where  coalescence 
is  complete,  scaling  using  combined  charge  weights  of  the 
coalesced  pulse  is  effective. 

Initial  peak  horizontal  velocity  similar  to  initial 
stress  is  the  result  of  input  from  charge  G  (20  tons)  at  all 
stations  except  the  deep  gages  300  and  400  feet  from  G. 

Peak  values  at  these  stations  scale  using  cube-root  scaling 
and  combined  charge  weights  of  F,  G,  and  H  (80  tons) .  These 
data  compare  well  to  Pre-GONDOLA  I  data. 

Peak  acceleration  predictions  for  the  crater  slopes  were 
based  on  competent  material  and  did  not  take  into  account  the 
presence  of  the  intervening  crater  rupture  zone.  This  method 
resulted  in  good  predictions  at  the  CP,  but  did  not  predict  the 
low  motions  on  the  Charlie  and  Delta  crater  slopes.  Maximum 
surface  accelerations  on  the  crater  slopes  were  strongly  influ¬ 
enced  by  the  effects  of  crater  boundaries  and  rubble. 
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APPENDIX  A 


STRESS-,  ACCELERATION-,  VELOCITY-, 
AND  DISPLACEMENT-TIME  HISTORIES 
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Figure  A.l  Stress  time  history,  Station  100 
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Figure  A. 2  Radial  velocity  time  histories  Station  100 
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Figure  A. 3  Stress  time  histories,  Station  150 
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Figure  A. 4  Vertical  velocity  time  histories,  Station  150 
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Figure  A. 5  Radial  velocity  time  histories,  Station  150 
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Figure  A. 6  Stress  time  histories,  Station  200 
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Figure  A. 7  Vertical  velocity  time  history,  Station  200 
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Figure  A. 8  Radial  velocity  time  histories,  Station  200 
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Figure  A. 9  Stress  time  histories,  Station  300 
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Figure  A. 11  Radial  velocity  time  histories,  Station  300 
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Figure  A. 12  Stress  time  history,  Station  400 
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Figure  A. 13  Radial  velocity  time  history,  Station  400 
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Figure  A. 14  Computed  radial  displacement  time  histories,  Station  100 
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Figure  A. 17  Computed  radial  displacement  time  histories,  Station  150 
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Figure  A. 18  Computed  vertical  displacement  time  history,  Station  200 
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Figure  A. 19  Computed  radial  displacement  time  histories,  Station  200 
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Figure  A. 21  Computed  radial  displacement  time  histories,  Station  300 
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Figure  A. 22  Computed  radial  displacement  time  history,  Station  400 
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